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Abstract. 

We present results from a comparative study of XMM- Newton observations of four classical T Tauri stars (CTTS), 
namely BP Tau, CR Cha, SU Aur and TW Hya. In these objects coronal, i.e. magnetic, activity and as recently 
shown, magnetically funneled accretion are the processes likely to be responsible for the generation of X-ray 
emission. Variable X-ray emission with luminosities in the order of 10 30 erg/s is observed for all targets. We 
investigate light curves as well as medium and high-resolution X-ray spectra to determine the plasma properties 
of the sample CTTS and to study the origin of their X-ray emission and its variability. The emission measure 
distributions and observed temperatures differ significantly and the targets are dominated either by plasma at 
high densities as produced by accretion shocks or by predominantly hotter plasma of coronal origin. Likewise the 
variability of the X-ray luminosity is found to be generated by both mechanisms. Cool plasma at high densities 
is found in all stars with detected O VII triplet emission, prevented only for SU Aur due to strong absorption. 
A general trend is present in the abundance pattern, with neon being at solar value or enhanced while oxygen, 
iron and most other metals are depleted, pointing to the presence of the inverse FIP effect in active coronae and 
possibly grain formation in evolved disks. We find that both accretion shocks and coronal activity contribute to 
the observed X-ray emission of the targets. While coronal activity is the dominant source of X-ray activity in the 
majority of the CTTS, the fraction for each process differs significantly between the individual objects. 
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1. Introduction 

Young, late-type pre main-sequence stars, known as 
T Tauri stars, are copiously found in or near star forming 
regions. Historically, they are classified according to their 
Ha equivalent width as classical T Tauri stars (CTTS) 
with EW>10A or otherwise as weak line T Tauri stars 
(WTTS). The principal physical difference between these 
two classes is that CTTS are thought to be in an earlier 
evolutionary stage, i.e. they still possess a disk contain- 
ing dust and gas and are accreting matter. WTTS are 
thought to be more evolved, to have mainly lost their disks 
and are approaching the main-sequence. While CTTS are 
generally younger with ages of a few up to ~ lOMyr and 
the overall fraction of CTTS in a given star-forming re- 
gion decreases with age, both types of T Tauri stars are 
commonly found in the same star-forming regions, indi- 
cating individual evolutionary time scales. Since the Ha- 
line is often time- variable, a classification solely based on 
its strength is sometimes misleading and other or addi- 
tional spectral features have been used to identify and 
classify CTTS more reliably. The existence of a disk con- 
taining a significant amount of matter and ongoing accre- 
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tion onto the host star impose major differences between 
the two types of TTS, leading to different spectral proper- 
ties which can be observed at different wavelengths. The 
effects of the warm and probably structured disk are for 
instance reflected in the different infrared designations of 
CTTS (Class II) and WTTS (Class III) due to their ad- 
ditional near-infrared emission, which is an important in- 
dicator for CTTS. Consequently, the different criteria and 
underlying data used for the identification of CTTS may 
introduce selection effect in the respective samples. 

Both types of T Tauri stars are well known 
strong and variable X-ray emitters and large num- 
bers of these objects were detected with E i nstein 
l|Feigelson fc DeCamplil. Il98lt iFeigelson fc Krisst Il989 ) 
and ROSAT l|Feigelson et all Il993t iNeuhauser et al , 
1199.1) . Their X-ray emission and the observed flaring was 
usually interpreted as a scaled-up version of coronal ac- 
tivity in analogy to cool main-sequence stars. This pic- 
ture is still valid and the X-ray emission from WTTS is 
thought to originate from coronal activity. However, T 
Tauri stars are thought to be at least in their early stages 
fully convective, therefore the question arises about the 
underlying processes for the generation of the magnetic 
field and the influence of accretion on the stellar struc- 
ture and corona. Further on, the evolutionary stages of 
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young stars and their high energy emission is of major 
importance for the process of planetary formation via its 
effects on chemistry, dust and protoplanets. A compre- 
hensive review of high energy processes and underlying 
physics in young stellar objects, summarising the results 
of the er a prior to Chandra and XMM-Newton, is pre- 
sented bv IFeigelson k Montmerld l|l999j) . 

The largest sample of pre main-sequence stars are 
found in the Orion Nebula C lo ud. U sing Chandra ACIS 
observations, iFeigelson et alJ l)2003l) analysed a sam- 
ple of 500 PMS stars with known basic properties. 
An even deeper exposure, named the Chandra Orion 
Ultr adeep Project (COU P) was presented recently, see 
e.g. lPreibisch et aP 1 2005|) . confirming and extending their 
results. The X-ray luminosity of CTTS in this sample 
was found to be correlated with bolometric luminos- 
ity. No strong correlation between Lx and rotation as 
observed for main-sequence stars was found, contrary, 
Lx correlates with Lboi and stellar mass, indicating the 
presence of different, e.g. turbulent, dynamos or of su- 
persaturation effects in solar-type dynamos. The ques- 
tion, whether the presence of an accretion disk influences 
the stellar X-r a y emi ssion is still under debate. While 
IFeigelson et al.l (j2008h found no intrinsic differences of 
X-ray activity related to th e presence of a circum stel- 
lar disk, the results found bv lPreibisch et all (|2005^ indi- 
cate that accretion diminishes Lx on average and weaken 
the Lx I Lh n ] correlation . A sim ilar influence was found 
by IStelzer fc Neuhauserl l)200l|) in the X-ray properties 
of young stars in the Taurus-Auriga complex observed 
with ROSAT PSPC. The COUP sources are rather con- 
sistent with the bulk of their X-ray emission being pro- 
duced in typical coronae. In general, magnetic activity in 
CTTS may also involve the circumstellar material, e.g. via 
star-disk interaction, but th e topology of their m agnetic 
fields in virtually unknown; iFavata et al.l l)2005|) present 
evidence for magnetic star-disk interaction in the analysis 
of larger flares on several COUP objects. 

Moreover, recent high-resolution spectroscopy of 
CTTS with Chandra and XMM-Newton, indicated ac- 
cretion shocks at least as an additional mechanism for 
the generation of X-ray emission in CTTS. In this in- 
terpretation X-ray emission is generated by a magnet- 
ically funnclcd accretion stream falling onto small ar- 
eas of t he stellar surface producing shocked high density 
plasma l|Shu et al l ll994tlCalvet k Cullbringl [l998). Since 
in magnetospheric accretion the matter falls from several 
stellar radii at nearly free-fall velocity, strong shocks and 
consequently hot plasma emitting in the X-ray regime are 
formed. This contribution is expected to be at energies 
below 0.5 keV and is distinguishable from cool coronal 
plasma by its density; therefore it can be recognised only 
through high resolution X-ray spectroscopy. It is thus not 
contradictory that the COUP sample based on ACIS data 
with medium resolution spectroscopy and poorer sensitiv- 
ity at low energies finds accretion to play no significant 
role. 



TW Hya was the first example of a CTTS, where obser- 
vational evidence for an accretion scenario could be found 
in X-ray data. Analysis using density sensitive lines of 
He-like triplets (e.g. O vn, Neix) indicated very high den- 
sities, exceeding by far the densitie s previously found in 
coronal plasma from any other star l|Kastner et all [2002; 
IStelzer k Schmittl . l2004j) . These lo w line ratios w e re als o 
found in the spectra of BP Tau l|Schmitt et all l2005|) . 
again strongly supporting an accretion scenario. However, 
accretion is only capable to produce plasma with temper- 
atures of a few MK and therefore contributes nearly ex- 
clusively to the low energy X-ray emission. The observed 
flaring and additional hard emission attributed to plasma 
with temperatures of several tens of MK in BP Tau and 
other CTTS require the presence of an additional X-ray 
generating mechanism such as magnetic activity. In this 
work we use the term 'coronal activity', but more complex 
phenomena like star-disk interactions cannot be ruled out. 
While TW Hya like objects where accretion is thought to 
be the dominant source of the X-ray emission are proba- 
bly rare, the BP Tau observation revealed that accretion 
as an additional X-ray generating mechanism may be more 
common. Another X-ray generating mechanism, shocks in 
strong Herbig-Haro outflows that may also produce a soft 
X-ray excess, was proposed for jet-driving CTTS (which 
resem ble Class I protostars) like DC Tau ( Giidel et all 
2005). We note that this model cannot explain the soft 
excess in the case of TW Hya or the other sample CTTS. 



We utilise XMM-Newton data from four X-ray bright 
CTTS to study the origin of their X-ray emission in ac- 
cretion shocks vs. coronal activity. We specifically use 
medium and high resolution X-ray spectra and spectrally 
resolved X-ray light curves to investigate the physical 
properties of our target stars BP Tau, CR Cha, SU Aur 
and TW Hya and to check the relevance of these mecha- 
nisms for the generation of the X-ray emission and vari- 
ability in these young stars. This work is the first com- 
parative X-ray study of CTTS using high and medium 
resolution spectra so far, therefore it is complementary to 
large sample studies of a specific star forming region such 
as COUP. Moreover, the data of CR Cha and partly of 
SU Aur and BP Tau is analyzed and presented here for the 
first time. Results derived from RGS spectra of BP Tau 
and the simultaneou s obtained UV data w ere presented in 
our previous letter ijSchmitt et alll2005h. while the da ta 
on TW Hya first presented bv lStelzer k Schmitd (|2004 is 
re-analysed for purposes of comparison in a manner iden- 
tical to the analysis of all other sources. The plan of our 
paper is as follows: In Sect.|21we describe the individual 
targets and previous X-ray results, in Sect. El the obser- 
vations and the methods used for data analysis and in 
Sect.0]we present and discuss the results subdivided into 
different physical topics, followed by a summary and our 
conclusions in Sect. El 
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2. The CTTS sample 

BP Tau is a CTTS with spectral type K5-K7 associ- 
ated with the Taurus- Auriga sta r forming region at a dis - 
tance of 140 pc as discussed bv IWichmann et alJ {l998). 
Optical and infrared measurements clearly show the typ- 
ical c haracteristics of excess emission produced by accre- 
tion ijGullbring et all Il998|) . BP Tau has a moderate in- 
clination angle of roughly 30° and its di sk was found to 
be quite compact (|Muzerolle et all 2003). X-ray emission 
from BP Tau with Lx ~ 10 30 erg /s was already m e asure d 
with the Einstein Observatory ^Walter fc KuhiL Il98ll) . 
the RASS luminosity was determined as 0.7 x 10 30 erg/s 
^Neuhauser et allll99^ . lGullbring et al.ll)l997|) presented 
an analysis of simultaneous optical and ROSAT observa- 
tions of BP Tau. Finding no correlation between optical 
and X-ray variability they concluded that the optical emis- 
sion is attributed to accretion while the X-ray emission 
arises from magnetically active regions. Initial results of 
our XMM-Newton observation of BP Tau with special fo- 
cus on the RGS data and density diagnostics with the 
Hc-l ike triplet of O vii w ere presented in our previous let- 
ter ijSchmitt et alll2005T) . where evidence for high density 
plasma and/or intense UV radiation is presented, support- 
ing the presence of an accretion funnel in BP Tau and a 
rough estimate of BP Tau's mass accretion rate as derived 
from simultaneous UV measurements was given. 

CR Cha with spectral type K2 is located in the 
Chamaeleon I cloud, a star forming region at a dis- 
tance of 140-150 pc with a measured R ASS luminosity of 
1.5 x 10 30 erg/s iFeigelson et allll993l). It was classified 



as CTTS via Ha, see e.g . iReirmrth et alJ l|199fifl . and ac- 
cording to lMeeus et al. I l)2003h it is dominated by small 
amorphous silicates, indicating large amounts of unpro- 
cessed dust. The archival XMM-Newton data of CR Cha 
is analysed and presented here for the first time. 

SU Aur is classified as G2 subgiant and also a mem- 
ber of the Taurus- Auriga star forming region, there- 
fore we adopt a distance of 140 pc. For this very lu- 
minous and variable object also a wider range of val- 
ues of its optical prope rties is given in literature, see 
e.g. iDeWarf et alJ (2003). It was sometimes classified as 
WTTS, but a clear IR-excess confirms its status as a 
CTTS, viewed nearly edge on, with a rather high accre- 
tion luminosity l(Muzerolle et all 12003') . thus causing the 
strongest absorption signatures in our sample. SU Aur is 
one of the brightest known CTTS, its strong X-ray emis- 
sion was already detected wit h the Einstein Observatory 
I Feigelson fc DeCa,mr>liLlT98T|) with an X-ray luminosity of 
3.0 x 10 30 erg/s. A simila r value (3.7 x 10 30 erg/s ) was de- 
rived from RASS data bv lNeuhauser et all (^995), making 
it the brightest target in our sample, especially when con- 
sidering the unabsorbed X-ray luminosity. Some results of 
the XMM-Newton observati on are presented i n an a naly- 
sis of star forming regions bv lPallavicini et all 1 2004h . who 
showed that SU Aur is both extremely hot and variable. 



TW Hya is classified as a K7-K8 star and is the near- 
est object in our sample. At a distance of only 56 pc it 
is located in the TW Hydra associatio n, a nearby young, 
but d iffuse stellar associations, see e.g. IZuckerman et alJ 
ll200ll) . TW Hya is viewed nearly pole on and due to its 
proximity and lack of obscuring clouds it can well stud- 
ied at all wavelengths. With an estimated age of nearly 
lOMyr it is al so one of the oldest known stars, still ac- 
creting matter. IKastner et alJ 1 19991) analysed ASCA and 
ROSAT data of TW Hya and found that a model with 
plasma temperatures of ~ 1.7 and ~ 9.7 MK and an X-ray 
luminosity of ~ 2 x 10 30 erg/s describes both observa- 
tions well. TW Hy a was also observed w ith the Chandra 
HETGS detector l|Kastner et all l2002h . providing high 
resolution spectra, which were modelled with an iron de- 
pleted and neon enhanced plasma and a emission measure 
distribution showing a sharp peak around 3 MK. The anal- 
ysis of density sensitive lines indicated plasma at high den- 
sities (log n e =13) and the inferred plasma properties were 
found to be consistent with the X-ray emission to be gen- 
erated by funneled mass accretion from the circumstellar 
disk. Several moderate 'flares' were detected on TW Hya 
during the ASCA and Chandra observation; however, the 
untypical shape of the flare light curve and the seemingly 
constant spectral properties led IKastner et al ] l|2002h to 
conclude that most or all of the X-ray emission of TW Hya 
originates from accretion. Results of the XMM-Newton 
observation of TW Hya obtain ed with a somewha t differ - 
ent analysis were presented bv lStelzer &: SchmittI l)2004|) . 
who found that the XMM-Newton data can be well ex- 
plained by X-ray emission from a metal depleted accre- 
tion shock with an X-ray luminosity of ~ 1.5 x 10 30 erg/s, 
consistent with the Chandra results. Different conclusions 
were derived at on the interdependent properties mass ac- 
cretion rate and surface area of the shock region, which 
is assumed to fill cither below one percent or up to a few 
percent of the stellar surface. Here the TW Hya data is 
independently re-analysed to ensure consistency and allow 
a comparison of our results. 



Table 1. Basic properties for our sample CTTS. 
Spectral types taken from Simbad da tabas e, dis- 
tances adopted from IWichmann et alJ lll998l). ages 
from iGullbring et"aD lll998h: B P Tau. iNatta et alJ 
ll20noh: CR Cha, IPeWarf et a,l J l|2003h : SU Aur. 

lMa.karov Rr. Fa.hricinsl l200ll) : TW Hva. L x are ROSAT 
values (< 2.4 keV). 



Target 


Spec.T. 


Dist. (pc) 


Age (yr) 




BP Tau 


K5V 


140 


6 X 10" 


0.7 


CR Cha 


K2 


145 


1 x 10 6 


1.5 


SU Aur 


G2III 


140 


4 x 10 6 


3.7 


TW Hya 


K8V 


56 


8 X 10 6 


2.0 
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3. Observation and data analysis 

All sample CTTS were observed with XMM-Newton using 
somewhat different detector setups with exposure times 
in the range of 30-130ks. Data were taken with all X-ray 
detectors, which were operated simultaneously onboard 
XMM-Newton, respectively the EPIC (European Photon 
Imaging Camera), consisting of the MOS and PN detec- 
tors and the RGS (Reflection Grating Spectrometer). Note 
that for SU Aur no PN data is available and different fil- 
ters were used, the thick filter for the BP Tau and SU Aur 
and the medium filter for the CR Cha and TW Hya ob- 
servations. Further on, the signal to noise ratio differs for 
the various targets and instruments. For CR Cha only a 
fraction of the original data from the additionally split ob- 
servation could be used for analysis because of very high 
background contamination; the RGS data from two ex- 
posures was merged with the tool 'rgscombine' provided 
with SAS 6.5, but the RGS spectrum of CR Cha is still 
underexposed. The EPIC data is quality is sufficient for 
all targets. In addition, strong absorption affects the data 
quality especially at lower energies. The details of our ob- 
servations and used data are described in Table |2 a de- 
tailed descr iption of the XM M-Newton instruments can 
be found in lEhle et all <|2003l) . 

Table 2. Observation log of our sample CTTS, duration 
of prime instrument/RGS (filtered). 



Target 


Obs.Modc 


Obs. Time 


Dur. (ks) 


BP Tau 
CR Cha 
SU Aur 
TW Hya 


FF/thick 
FF,LW/mcd. 
FF/thick 
FF/mcd. 


2004-08-15T06:14-16T18:51 
2001-02-24T05:04-15:49 

2001-09-21T01:27-22T14:18 
2001-07-09T05:51-14:01 


132/124 
39/72 

130/126 
30/28 



Data analysis was performed with the XMM-Newton 
Science Analysis System (SAS) software, version 6.0. 
Images, light curves and spectra were produced with stan- 
dard SAS tools and standard s election criteria were ap- 
plied for filtering the data, see lLoiseau et all ((2004) . X- 
ray spectral an alysis was carried out with XSPEC VI 1.3 
ijArnaudl Il996|) . whil e for line fitting pu r poses we used 
the CORA program l|Ness fc WichmannL 12002^ . assum- 
ing Lorentzian line shapes. Individual line fits are used 
to investigate the density sensitive forbidden and inter- 
combination lines of He- like triplets. Line counts are de- 
rived keeping the line spacing fixed within a triplet and 
using an overall line width for all lines. Spectral analy- 
sis of EPIC data is performed in the energy band between 
0.3- 10.0 keV, while the RGS first order spectra in the full 
energy range, i.e. 0.35-2. 5keV (5-35A), are used when- 
ever data quality permits. While the RGS obviously has 
the highest spectral resolution, the EPIC detectors are 
more sensitive and are able to measure higher energy X- 
rays; the MOS detectors provide a slightly better spectral 
and spatial resolution, the PN detector is more sensitive. 
We emphasize that the data were analysed simultaneously 
but not co-added, thus providing sufficient signal for all 
observations and ensuring a consistent analysis for all tar- 



gets. All periods affected by proton flares were removed 
from spectral analysis and for TW Hya the PN data was 
cleaned for some moderate pile-up. The background was 
taken from source free regions on the detectors. 

For the analysis of the X-ray spectra we use multi- 
temperature models with variable but tied abundances, 
i.e. the abundance pattern was assumed to be the same 
in all temperature components. Such models assume 
the emission spectrum of a collisionally-ionized optically- 
thin gas as calcu lated with the APEC code, see e.g. 
Smit h et al ] ll200lh . Since the RGS is not very sensitive 
above 2.0 keV, we use EPIC data to determine the prop- 
erties of plasma with temperatures above 20 MK. In order 
to account for calibration uncertain ties of the different de- 
tectors, see e.g. lKirsch et alJ l|2004l) . the normalization be- 
tween instruments was taken as a free parameter for each 
type of instrument, i.e. for MOS, PN, RGS. Our fit pro- 
cedure is based on \ 2 minimization, therefore spectra are 
always rebinned to satisfy the statistical demand of a min- 
imum value of 15 counts per spectral bin. All errors are 
statistical errors given by their 90% confidence range and 
were calculated separately for abundances and tempera- 
tures by allowing variations of normalizations and respec- 
tive model parameters. Note that additional uncertainties 
arise from uncertainties in the atomic data and instrumen- 
tal calibration which are not explicitly accounted for. 

The applied models use three temperature compo- 
nents, models with additional temperature components 
were checked, but did not improve the fit results signif- 
icantly. Since the APEC models do not account for high 
density plasma, resolved density sensitive lines are ex- 
cluded from the global fits for TW Hya. Abundances are 
calculated relative to solar photospheric values as given by 
lAnders fc Grevessel lll989l). For iron and oxyg en we use the 
updated values of iGrevesse fc SauvaH l|l998lh Application 
of the new but co ntroversial solar abundances published 



hv lAsnhmdl 420051) would further increase the neon abun- 
dance in our sample stars compared to other metals, but 
not change our conclusions significantly. We note that rel- 
ative results between the CTTS are not affected by the un- 
derlying solar abundance pattern. For elements with over- 
all low abundances and no significant features in the X-ray 
spectra, i.e. Al, Ca, Ni, the abundances were tied to the 
iron abundance. When data quality permits, we determine 
the abundances of individual elements, global abundances 
are used if no features are present in the spectra. This is 
for instance the case for targets with strong absorption, 
where e.g. the stronger carbon lines at the low energy end 
of the detectors are completely absorbed. We simultane- 
ously modelled the temperatures and emission measures 
(EM= J n e nffdV) of the components and checked the de- 
rived results for consistency. X-ray luminosities were then 
calculated from the resulting best fit models. 

Absorption in the circumstcllar environment and pos- 
sibly also in the interstellar medium is significant for 
CTTS and is applied in our modelling. Since absorption is 
supposed to be slightly variable and, in addition, optical 
and IR measurements often give a wider range of values, it 
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is kept as a free parameter in our analysis. In general, the 
derived fit results are quite stable, but note that some of 
the fit parameters are mutually dependent, thus affecting 
especially absolute values. Interdependence mainly affects 
the low energy region of the spectra, where the strength 
of absorption, the emission measure and abundances of el- 
ements with emission lines in the respective temperature 
range strongly depend on each other. While some effects 
can be prevented by the methods described above, models 
with different absolute values of the mentioned parameters 
but only marginal differences in its statistical quality may 
be fitted to the data, however ratios and relative changes 
of these properties are again very robust. 

4. Results 

4.1. X-ray light curves 

To study time variability and its origin we first investi- 
gated the X-ray light curves of our sample CTTS. In Fig|T] 
we show background subtracted light curves with a tempo- 
ral binning of 1000s, extracted from a 50 "radius circular 
region around each source and cleaned for obvious data 
dropouts. 

A flare with a rise in count rate of factor 2.5-3 oc- 
curred during the BP Tau observation and a smaller 
flare (factor ~1.5) during the CR Cha observation. 
In the SU Aur data several flares were detected, the 
largest one showing an increase in X-ray brightness by 
a factor of 3. TW Hya does not exhibit a flare dur- 
ing the XMM-Newton observation, any variability re- 
mains in the range of only 10%. However, light curve 
variations with factors around two a re also known for 
TW Hya l|Kastner et all Il999t l2002j) . who noted three 
'flares' during 140 ks (ASCA, 94 ks + Chandra, 48 ks) ob- 
servation time but no spectral changes were detected dur- 
ing the periods of increased X-ray brightness. The average 
time scale for flares on CTTS appears to be around one 
moderate flare (factor > 2) occurring daily on our sample 
CTTS. In addition to obvious flare events, variability is 
present also on smaller time scales and lower amplitude 
variations throughout the observations. 

4.1.1. X-ray hardness 

The light curves in connection with spectral hardness can 
be used to identify the origin of the variability. While 
in typical stellar coronal flares the emission measure of 
predominantly hot plasma and hence the hardness of the 
spectra is increased, in a pure accretion spectrum no spec- 
tral changes should accompany the brightening, since the 
plasma temperature only depends on the infall velocity 
and not the accretion rate. If, in addition, a coronal con- 
tribution with a temperature higher than that produced 
by accretion is present, a slight spectral softening should 
be observed. We calculate a hardness ratio for each light 
curve time bin and in Fig-EI w ^ show the hardness ratio 
vs. count rate for our sample CTTS. The hardness ratio is 



i BP Tau 
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Time (ksec) 



: CR Cha 




0.1 5 L 



10 20 30 
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SU Aur 




0.2 - 



20 40 60 80 100 120 

Time (ksec) 



TW Hya 




2.-1 h - 

5 10 15 20 25 

Time (ksec) 

Fig. 1. Light curves and corresponding hardness ratio of 
the sample CTTS, PN (SU Aur-MOS) data with lksec 
binning. 

here defined as HR=H-S/H+S with the soft band cover- 
ing the energy range 0.2-1.0keV and the hard band 1.0- 
10.0 keV. Errors are small compared to the observed shifts 
in hardness ratio. A clear correlation of X-ray brightness 
with spectral hardness is present for BP Tau, CR Cha and 
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SU Aur, a behaviour typical of stellar flares, suggesting 
a coronal origin of the variability. No strong correlation 
is found for variations on TW Hya, which are, however, 
quite small. An anti-correlation appears to be present for 
the larger variations, as expected for brightness changes 
due to an increase in accretion rate. This is reflected in 
the slope of the linear regression curves. It is positive for 
BP Tau, CR Cha and SU Aur, but slightly negative for 
TW Hya. 
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Fig. 2. Hardness ratio vs. count rate of the sample CTTS, 
PN (SU Aur-MOS) data with lksec binning. 



To quantify the spectral changes that the CTTS are 
undergoing during phases of increasing emission in more 
detail, we use their light curves to define 'high' states for 
BP Tau (15-37ks), CR Cha (18.5-31.5ks), SU Aur (0-24, 
47-69, 83-101 ks) and TW Hya (> 19 ks), while the re- 
maining time intervals are considered as quasi-quiescence. 
An individual spectral analysis of both states of activity 
is presented in Sect. 14.2.31 



4.2. Spectral analysis 

The spectral analysis and its results are presented as fol- 
lows. First we use the medium resolution EPIC spectra 
in combination with the high resolution RGS data to ob- 
tain general spectral properties like temperatures, emis- 
sion measures and elemental abundances of the X-ray 
emitting plasma and constrain the strength of absorption 
present for our targets. From the RGS spectra we investi- 
gate the resolved He-like triplet of O vn to determine the 
densities of the emitting plasma and/or strength of the 
surrounding UV field. Subsequently we investigate spec- 
tral changes between the low and the high state of our 
sample CTTS. 



4.2.1. The global spectra 

The spectra of all our sample CTTS as observed with 
MOS1 are shown in Fig.|3 which demonstrates major dif- 
ferences between individual stars. While the spectra of 
BP Tau and CR Cha are comparable, the spectrum of 
SU Aur indicates the presence of large amounts of ex- 
tremely hot plasma noticeable, e.g. in the very strong 
FeXXV line complex at 6.7 keV. On the other hand, the X- 
ray spectrum of TW Hya is much softer, suggesting a more 
dominant accretion component. The observed spectra are 
also subject to absorption, affecting primarily the energy 
range below 1.0 keV. In these spectra the low energy slope 
mainly reflects the strength of the absorption, while the 
high energy slope traces temperature and amount of hot 
coronal and flaring plasma. Inspection of the two slopes 
indicates that the absorption is weakest for TW Hya, mod- 
erate for BP Tau and strongest for CR Cha and SU Aur, 
while the coronal component is strongest and hottest for 
SU Aur, followed by BP Tau and CR Cha and much 
weaker and cooler for TW Hya. 




Energy (keV) 

Fig. 3. Spectra of our sample CTTS for the total observa- 
tion as observed with the MOS1 detector: BP Tau (black), 
CR Cha (green), SU Aur (blue) and TW Hya (red); bot- 
tom to top of the 1.0 keV peak. Colour figure in electronic 
version. 

To quantify the physical properties of the X-ray emit- 
ting plasma we fitted the spectra of the different in- 
struments with multi-temperature spectral models as de- 
scribed in Sect |21 The results of the fitting procedure are 
presented in Tabled! Note again the interdependence of 
several fitting parameter, especially absorption, cool emis- 
sion measure and corresponding abundances. However, 
general trends on abundances and their ratios as well as 
the overall shape of the emission measure distribution and 
their changes, which indicate the production mechanism 
of the X-rays, are very stable. 

While spectral modelling must consider absorption, 
the physically most interesting quantity is the emitted, 
i.e. dereddened, X-ray luminosity of our targets, which is 
calculated in the 0.2-10. keV band using MOS data (see 
01) • We find that despite their very different spectral prop- 
erties, all targets have an X-ray luminosity of the same 



Robrade & Schmitt: CTTS in X-rays 



7 



order of magnitude. The hottest and most massive star 
SU Aur is also X-ray brightest target, followed by CR Cha, 
BP Tau and TW Hya with comparable X-ray luminosities. 
CTTS as a class are known strong X-ray emitter, but our 
sample stars are about a magnitude brighter than the X- 
ray luminosity of an average CTTS, as determined from 
RASS data with a mean logLy of 29.1 (Chamaeleon) and 
29.4 ( Taurus) in iNeuhauser et alJ l)l995|) . lEsigelson et"aT1 
( 1993J); this is of course no surprise since these CTTS were 
selected for grating observations. 

Table 3. Spectral results, units are Nh in 10 21 cm -2 , kT 



inkeV EM in 10 5 



in 10 ergs 



n h 

Fe 

Si 

O 

Ne 

kTl 

kT2 

kT3 

EMI 

EM2 

EM3 



0.287 
0.14+ 



BP Tau 

1 5+0-1 

i -°-o.i 

3 + 0.04 
-0.04 
4+0.08 
-0.0S 
n co + 0.08 

1 47+ - 16 

-0.14 

o 20+ 001 
u -^ u -o.oi 

fi3+ 03 
UDJ -0.03 
9 , 7 +0.09 

3 48+ 89 

°-*°-0.69 

5 28+ 0134 
D - ZS -0.56 
-, + 0.43 
-0.29 



10.30^ 



CR Cha 

9 7 + 0-4 

z - ' -0.3 

45+ - 08 
23+ 10 

R3+ - 18 
u -°°— 0.14 

1 02 +0 19 
17+° 03 

u - 1 ' -0.03 

0.66I 02 

+ 0.06 
0.10 
+4.05 
2.49 
+ 1.60 
1.53 
+0.73 
0.54 



1.92 
4.96 
10.67 
9.27 



SU Aur 
3 n+0.1 — 
O,i -o.i 
n 7q+o.05 

u - ' "-0.05 

0-60i o°o 8 8 

1 04 +01B 
n fis +o.03 

°- faf '-0.02 
, +0.07 
L -0.10 
-,+ 1.02 
-0.84 
- + 0.87 



1.61] 
4.82 4 



H-66-0.86 
25.27l 2 3 ;^ 
15 79+ 4 02 

LO - ' -2.59 



TW Hya 

u - oo -0.05 
n nc + 0.02 
°- 26 -0.02 
n nn + 0.05 
u - 1J -0.06 
„ „„ + 0.02 
u - ZJ -0.02 
, 01 +0.05 
!- 81 -0.05 
- + 0.01 
-0.01 

69+ 06 

u - Da -0.02 

1 34+ 03 
19 16+ 074 

, ,,+0.54 
L -* o -0.8i 
, + 0.52 
-0.51 



0.25! 



3.12"* 



X A (d.o.i.) 


0.97 (1281) 


1.20 (700) 


1.20 (1027) 


1.57 (878) 


Lx obs. 


1.3 


1.4 


4.5 


1.5 


Lx emit. 


2.3 


3.1 


8.1 


2.0 



Neglecting interstellar contribution, the strength of 
absorption depends on the amount of circumstellar ma- 
terial in the disk and the inclination angle of the tar- 
get. The results of our spectral modelling confirm that 
absorption is lowest for TW Hya. Our value for the in- 
terstellar colu mn density Nh is s lightly higher than the 
one used by IStelzer fc Srlmiinl 0004), who fixed N H 



at 2 x 10 u cm . No values are given for the Chandra 
data in iKastner et al.| ||2002|). but the authors refer to 
previous work |KasTrieretTll 1999) where values deter- 
mined from ROSAT PSPC (0.5 x 10 
ones from ASCA SIS (2.9 x 10 



21 



21 



cm 
- 2 ) 



~ 2 ) and higher 
with a conse- 
quently larger emission measure. These differences were 
attributed to temporal changes in the circumstellar envi- 
ronment, but the poorer spectral resolution and sensitiv- 
ity of the used instruments, the strong interdependence 
between Nh and EM, that is also present in models from 
the individual instruments operated simultaneously on- 
board XMM-Newton, suggests that limited spectral reso- 
lution of the instruments combined with uncertainties in 
the code used for modelling the spectra is a possible ex- 
planat ion for those discrepancies. Recently iHerczeg et alJ 
( 2004) analysed the Ly a profile of TW Hya as measured 
by HST/STIS and put a stringent upper limit on the hy- 
drogen column density with Nh < 5 x 10 19 cmT 2 . We 
therefore adopted this value and additionally fitted the 
spectra with the reduced absorption value, which results 
in a comparable model, but with accordingly lower emis- 
sion measure at cool temperatures. 



We find intermediate values of the interstellar col- 
umn density Nh for BP Tau, which are in good agree- 
ment with optical measurements using the standard con- 
version between Nh and A v viz. Nh=2 x 10 21 A V cm~ 2 ; 
values take n from literature are in the range 1.0 — 2.0 x 
10 21 cm- 2 dOullbring et all Il998t Eartigan et all Il995h . 
Even stronger absorption is found for CR Cha, also in 
good agreement with optical me asurements, with aliter- 
ature value of 2.7 x 10 21 cm" 2 (jOauvin fc StromI Il992h . 
The strongest absorption is found for SU Aur, here the 
large inclination angle results in a strong and probably 
also variable absorption. Optical and NIR measurements 
give a range of values ranging from no absorption up to 



ot values ranging trom no absorption up to 
2 l|DeWarf et all 1200.^ . While commonly a 



3 x 10 21 . 

lower Nh value is adopted, our result is comparable with 
the larger values and a much weaker absorption signifi- 
cantly worsens the quality of the fit. 

The emission measure distributions (EMD) of our sam- 
ple CTTS differ significantly among the sample stars. The 
hottest EMD is clearly present in SU Aur, with significant 
amounts of plasma at temperatures around 50 MK de- 
tected. Due to the very strong absorption no definite con- 
clusions about a possibly existing cool temperature part of 
the EMD can be drawn. Since SU Aur is of average age in 
our sample and the major difference to the other stars is its 
spectral type, an explanation for its outstanding proper- 
ties might be a dependence of the evolutionary time scale 
on spectral type, i.e stellar mass. The EMDs of BP Tau 
and CR Cha are similar, with the temperatures of BP Tau 
being slightly higher and a larger fraction of plasma re- 
siding in the hottest component; plasma at temperatures 
around 20-25MK must be present in those stars. The 
picture is very different for TW Hya, where the cool com- 
ponent with temperatures around 3 MK clearly dominates 
the EMD. Plasma at medium temperatures with 5- 10 MK 
appears to be nearly absent, while in the hot component 
plasma temperatures around 15 MK are reached; however, 
this component is cooler and significantly weaker than in 
the other objects. We emphasise that this finding is inde- 
pendent of the used value of Nh- 

As far as abundances are concerned, one property is 
common to all CTTS analysed: Most abundances are sub- 
solar, sometimes at considerable level, while neon is com- 
monly found at solar abundance or even significantly en- 
hanced. The noble gas neon is more abundant compared 
to iron, oxygen and silicon for all targets, and in the 
most extreme case TW Hya these ratios are of the or- 
der of 10 (see Table. 01 • Iron and silicon are less abundant 
than oxygen in BP Tau and CR Cha, while all three el- 
ements are strongly depleted in TW Hya. The very ac- 
tive star SU Aur appears to be an exception, with only 
oxygen apparently being more strongly depleted, but here 
very strong absorption affects the results based on cool 
lin es. These metal anomalies were interpreted for TW Hya 
bv lStelzer k. Schmittl l)2004|) through a depletion of grain 
forming element in the accreted material, consistent with 
the low observed IR excess an d also plausible wh en con- 
sidering the age of TW Hya. iDrake et alJ l|2005h argue, 
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that for TW Hya coagulation of grains into larger bod- 
ies finally withdraws these elements from the accretion 
process, while for the younger BP Tau dust and grains 
at corotation radius (R co = (^r-) 1 / 3 ) will sublimate and 
rereleased into the accretion process. The coagulation of 
significant amounts of material into centimeter sized par- 
ticles in the disk around TW Hya is supported by the 
radiospec trum at centimeter w avelength as observed with 
the VLA iWilner et al 112(1051) . 

However, stars without accretion but with an active 
corona also show a distinct abundance pattern, com- 
monly known as IFIP effect, i.e. an enhancement of el- 
ements with a high first ionisation potential. While iron 
and silicon are low FIP elements, oxygen is an inter- 
mediate and neon a high FIP element. The IFIP-effect 
is observed for BP Tau and CR Cha, thus pointing to- 
wards the typical abunda n ce pa ttern for coronal plasma. 
Moreover, lArgirofii et alJ l|2005|) analyse the spectra of 
TWA5, a young multiple system dominated by a M dwarf 
and also located in the TW Hya association. These ob- 
ject is classified as CTTS (via Ha), but is probably a 
WTTS since there is no evidence for a disk in the IR, but 
it also shows a high neon to iron ratio similar to TW Hya 
with spectral properties very much reminiscent of active 
stars. While environmental conditio ns may be invoked 
to exp lain the abundance anomalies, iR.obrade fc Sch mitt 
;'2( )(),") • showed that a high neon to iron ratio is common for 
active M dwarfs, with most extreme ratios found in very 
young and therefore more active stars. We further point 
out that our CTTS belong to different star formation re- 
gions and their X-ray emission is predominantly generated 
by different processes, yet no significant difference in the 
abundance pattern was found. Therefore it appears rea- 
sonable to attribute the observed abundance anomalies in 
the sample CTTS to a combined action of metal depletion 
of the accreted matter via grain forming and the coronal 
IFIP effect. Grain forming is only important for the ac- 
cretion plasma of older and more evolved objects, while 
the IFIP effect is present in plasma generated in active 
coronae. 



4.2.2. The Ovm and Ovu lines 

We used the CORA line fitting program to determine the 
strengths of the resonance, intercombination and forbid- 
den lines in the He-like triplet of O vn (21.6, 21.8, 22.1 A) 
and the Lya line of Ovm at 18.97 A as measured with 
RGS1. The spectra of the sample CTTS covering this 
wavelength-region are shown in Fig.^ the measured line 
counts of the O vn triplet and the derived densities are 
given in Tabled We used the relation f/i= i+^/^Xne/N 
with f and i being the line intensities in the forbidden and 
intercombination line, Rq the low density limit of the line 
ratio with a adopted value of 3.95, N c the critical density 
and 4>/ 4> c the radiation term, which is neglected in our cal- 
culations. Values used in the calculations were taken from 
I Pradhan & Shulll Il98llk we caution that the presence 



of strong radiation fields would lead us to overestimate 
plasma densities. The O vn triplet traces the cool plasma 
around 2MK, is essentially free of stronger blends and 
therefore well suited to investigate possible accretion sce- 
narios, while O vni traces slightly hotter plasma around 
3-5MK. 

The analysis of the oxyg en triplet of BP Tau and its 
implications are presented bv lSchmitt et al J (120051). results 
obtain ed for TW Hya are presented in lStelzer fc Schmittl 
(2004). Our analysis here is - technically - slightly dif- 
ferent, but we arrive a t the same results as found by 
IStelzer fc Schmittl l|2004|) . The data of CR Cha, which un- 
fortunately has only moderate S/N ratio and larger ab- 
sorption, is analysed here for the first time. The O vni 
line is here clearly detected, the same applies to the O VII 
triplet, which, however, looks quite peculiar. Resonance 
and forbidden lines have about the same strength, and 
the intercombination line is actually the strongest line in 
the triplet. Because of the poor SNR the errors in the 
line count measurements are quite substantial; also the g- 
ratio (f+i)/r is found to be 2.9 ±2.2, but is because of its 
large error actually consistent with unity as theoretically 
expected. The observed f/i ratio is well below unity de- 
spite its considerable measurement error, as observed for 
TW Hya and BP Tau. 

In SU Aur the Ovu triplet is absent, ruling out any 
quantitative analysis. We attribute the absence of Ovu 
to the strong (circumstellar) absorption and the pres- 
ence of extremely hot plasma, which increases the con- 
tinuum emission at these wavelengths. Some weak emis- 
sion may actually be present at the position of the oxy- 
gen triplet, but neither the intercombination nor forbidden 
line are clearly detected. The Ovm Lya line is however 
clearly present and attenuation increases only by a fac- 
tor of roughly two between the two wavelengths. A strong 
Ovu triplet should therefore be detected in the spectra, 
even with SU Aur's strong absorption. This indicates, that 
large amounts of cool plasma do not exist on SU Aur. 

Table 4. Measured RGS1 line counts of the Ovm and 
Ovn-triplet (r,i,f) lines and calculated densities, taking 
into account the effective areas at the respective wave- 
lengths. 



Par. 


BP Tau 


CR Cha 


SU Aur 


TW Hya 


O vin 


94.5±11.4 


36.6±7.1 


33.1±7.5 


453.0±23.8 


O vn r 


44.9±8.4 


7.5±4.9 




234.9±17.3 


O VII i 


34.5±7.6 


13.5±5.6 




113.5±12.9 


O vn f 


11.6±5.3 


8.7±4.8 




6.1±5.1 


f/i (obs.) 


0.34±0.17 


0.64±0.44 




0.054±0.045 


n e (W xl cmT' 6 ) 




'■■"-0.8 







Comparing the ratio of O VIII to the O VII line counts, 
we find the largest (observed) ratio for TW Hya. However, 
absorption attenuates the O VII lines more strongly than 
the O vni line and using our determined values for Nh , in 
the unabsorbed ratio the inverse trend is observed. While 
typical active coronae, as well as all other sample CTTS, 
show a rise of the EMD in this temperature regime, we 
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Fig. 4. Spectra of the sample CTTS covering the Vin 
line and the Ovu triplet as measured with the RGS1 de- 
tector. Absorption increases from the top to bottom. 

find that in the accretion dominated spectrum of TW Hya 
the bulk of plasma has temperatures around 2.5- 3 MK, 
consistent with the results of the global fitting. 

Investigating the O vn triplets of our sample CTTS, we 
find that derived densities differ from those of cool main 
sequence stars. While typical coronal plasma is compati- 
ble with and for very active stars no more than an order 
of magnitude above the low den sity limit, i.e f/i ~ 2-4 and 
log n e ~ 9-10 l|Ness et all .20021) . all CTTS with detected 
oxygen triplets deviate strongly from the low density limit. 
BP Tau and CR Cha show a density about two orders of 
magnitude lower, for TW Hya it is even three orders of 
magnitude lower, suggesting extremely high densities for 
the cool plasma of the CTTS. Alternatively, strong UV 
radiation fields would explain the measured line ratios, 
but these fields cannot be produced by the stellar photo- 
spheres of the underlying cool stars, therefore the radia- 
tion fields would have to be attributed to an additional 
hotter component, i.e. an accretion shock. 



Thus the derived densities may be considered as an 
upper limit, but we point out that even in the case of 
UV radiation contributing to the low f/i-ratio accretion 
must occur, since also a strong UV radiation points to the 
presence of an accretion-induced shocks. Furthermore, the 
coronal plasma is likely to additionally co ntribute to the 
flux in the O vn-triplet. It was shown in ISchmitt et alj 
(2005J) that the f/i-ratio for BP Tau is even lower, when 
only the quasi-quiescent state is considered in the analy- 
sis. This is plausible, because the f/i-ratio from the coronal 
contribution is expected to be near the low density limit 
and therefore the densities in the accretion spot are under- 
estimated in this calculation and also appear as an lower 
limit. Since the size of the accretion spots or filling factors 
as well as the amount of coronal plasma contributing at 
low temperatures are not known precisely, we just note 
that there are two competing effects on the derived den- 
sities present, but they are not calling the accretion sce- 
nario in question. While the poor data quality for CR Cha 
prevents a more detailed discussion as for BP Tau, the de- 
rived f/i-ratio is also significantly lower than for coronal 
sources and the derived plasma density is in the same or- 
der of magnitude as it is for BP Tau. The derived results 
for CR Cha fit into the picture and support the accretion 
scenario for CTTS as a class. Therefore we conclude that 
the cool plasma of TW Hya is strongly dominated by ac- 
cretion, while in the cases of BP Tau and CR Cha and 
possibly also of SU Aur accretion is present, but with ad- 
ditional cool coronal plasma. However, the emission mea- 
sure of this component is hard to constrain given the SNR 
of our data. These findings agree well with the conclusions 
drawn in the previous section. 

4.2.3. Spectral changes 

To investigate the spectral changes between the phases 
of different X-ray brightness, we use the more sensitive 
EPIC data, separated in low and high states as defined in 
Sect l4.1l and utilise the previously derived best fit models. 
A comparison of the PN spectra of BP Tau and TW Hya 
during high and low state is shown in Fig.El where the 
high state represents the strong flare on BP Tau (fac- 
tor 2.4) and the period of slightly enhanced emission of 
TW Hya (factor 1.1). The differences between the spectra 
of the high and low state are nearly reverse for the two 
CTTS. While for BP Tau the spectrum hardens and the 
most significant changes are seen at higher energies, point- 
ing to a coronal origin of the changes, the TW Hya spectra 
are identical at higher energies and only at lower energies 
the emission is enhanced, pointing to a higher accretion 
rate. The changes are more pronounced for BP Tau due 
to the much higher increase in X-ray brightness. 

No significant changes in the elemental abundances 
were detected between high and low state, therefore we 
fixed strength of absorption and abundances at those val- 
ues that were derived for the total observation and remod- 
elled the EMD for the selected time intervals. If we find - 
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Fig. 5. Spectra of BP Tau (top) and TW Hya (bottom) 
during high (black) and low (red) state with model fits. 

Table 5. EMD of our sample CTTS in low (left) and high 
(right) state as derived from PN data (MOS for SU Aur), 
Njj in 10 21 cm- 2 , kT inkeV and EM in 10 52 cm- 3 . 



Flaring affects also the emission measure at moderate tem- 
peratures as represented by the increased emission mea- 
sure and a shift to higher temperatures, which are in the 
range of 6-8 MK. A more moderate increase of temper- 
atures and additional plasma at moderate temperatures 
but again predominantly in the hot component is found 
for CR Cha. This also points to a coronal origin of the 
increased brightness and is consistent with the observed 
more moderate flare. Also for SU Aur coronal activity is 
the explanation for the increased X-ray brightness, with 
large amounts of additional hot plasma at temperatures of 
20-50MK. The flare plasma is extremely hot, but since 
large amounts of extremely hot plasma are already present 
during quasi-quiescence, no significant increase in temper- 
ature accompanies the brightening. The coolest detected 
component of SU Aur, which already has temperatures 
around 8 MK, is again not strongly affected by the flaring. 
In contrast, none of these signatures are connected with 
the flux increase for TW Hya. No differences in the tem- 
perature structure were found and the additional plasma 
is found mainly in the coolest temperature component, 
supporting the scenario of increased accretion. With 2.5 — 
3MK this component is much cooler than typical flare 
plasma temperatures, but fits to the temperatures as ex- 
pected for plasma generated in an accretion spot. 



5. Summary and Conclusions 



Par. 



BP Tau 



CR Cha 



N H 

kTl 

kT2 

kT3 

EMI 

EM 2 

EM3 



1.5 



0.18 

u - ou -0.05 
r, -I , + 0.15 
— 0.14 

4.5 ' 



5.9 
6.8 



+0.5 
0.5 
0.4 
-0.5 
0.5 
-0.4 



+ 0.02 
0.02 
n c >7+0.07 

+ 0.26 
0.14 
+ 1.0 
-1.0 
+0.9 
-1.2 



2.71] 
2.! 
9.8"* 



2:;. 9 



+ 1.4 
-1.7 



1.71 
4.5 

10 
7.7 



2.6 

19+ 006 

u ' it '-0.04 
+ 0.03 
0.03 

1.99 



0.65 

+ 0.20 



0.19 
+ 1.5 
-1.2 
. + 1.5 
'-1.7 
+ 0.9 
-0.9 



+ 0.23 
0.14 
+ 2.9 
2.3 

11 7+ 1 ' 4 

n, '-1.9 
+ 13 
1.4 



6.3: 



13. 3 1 



X'(dof) 


0.89 (778) 


1.12 (414) 


Par. 


SU Aur 


TW Hya 



N H 

kTl 

kT2 

kT3 

EMI 

EM2 

EM3 

x"(dof) 



0.69 



3.1 

+ 0.03 
0.02 



10.8 
21.2 
6.9 



1.627 
5.08 

-1.0 
-0.5 
+ 1.8 



+ 1.9 
1.9 



+ 1.27 
-0.81 

11.9 
29.9 
22.7 



+ 1.1 

0.3 
+ 4.1 
4.3 
+ 4.4 
4.3 



0.35 

24+ 01 
u - z -o.oi 

60+ 09 

U.DU_q 06 

, oy + O.ll 

101 -0.05 

21.51°-* 24.51 

„ q + 0.8 o 
z -°-0.4 °'^-0.3 
3 + 0.3 r, O+0.5 



-+0.4 



-0.6 



2.3" 



1.10 (929) 



1.27 (496) 



using temperature as a free parameter for each state - no 
significant difference in the temperature components, it is 
tied for both sets of data. 

The results of the individual modelling of the low 
and high states for our sample CTTS are summarised in 
Tabled Very different effects on the derived temperature 
structures are found and also the largest increase in emis- 
sion measure is found in different temperature components 
for our CTTS. The rise in temperature and the presence 
of a larger amount of additional hot plasma with tempera- 
tures around 30 MK during the high state reflects the coro- 
nal origin of the increased X-ray brightness for BP Tau. 



We have presented the first comparative study using high 
and medium resolution X-ray spectra of classical T Tauri 
stars observed with the new generation X-ray telescopes 
so far. The results derived from these data are complemen- 
tary to the ones obtained from the large sample of COUP 
sources, since high-resolution spectroscopy and low-energy 
sensitivity additionally permits the investigation of the 
cool plasma components. 

Using the XMM-Newton observations of BP Tau, 
CR Cha, SU Aur and TW Hya we determined X-ray prop- 
erties of accreting young pre-main sequence stars. Two dif- 
ferent mechanisms are likely to contribute to the produc- 
tion of X-ray emission in those objects, first coronal, i.e. 
magnetic activity, and second, magnetically funneled ac- 
cretion. We investigate variability, global spectra and den- 
sity sensitive lines to check for the presence and relative 
contribution of these two mechanisms for the individual 
objects. For this purpose we utilize the emission measure 
distribution and its changes, specific abundance pattern 
and density diagnostics. Interpreting very high densities 
as indicator for accretion shocks and a high temperature 
component as indicator for coronal activity, we derive the 
following conclusions which are likewise supported by the 
analysis of abundance analysis and spectral variability. 
In all targets where the Ovil-triplet is observable den- 
sity analysis leads to higher densities than observed in 
any pure coronal source. Additionally, a high temperature 
component is present in all targets. We therefore argue 
that both X-ray generating processes are present in our 
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sample CTTS, but at very different levels and importance 
in individual objects. 

The X-ray emission of BP Tau is overall dominated 
by coronal activity and a stronger flare of again prob- 
ably coronal origin is observed. Even in quasi-quiescent 
phases the emission measure distribution is dominated 
by medium and hot temperature plasma whereof signifi- 
cant amounts are present at temperatures around 20 MK. 
However, the analysis of density sensitive lines confirms 
that accretion is also present and actually dominates the 
cool plasma with temperatures of a few MK. CR Cha is 
similar to BP Tau but more strongly absorbed and the 
coronal plasma temperatures are slightly lower. It also ap- 
pears to be dominated by the coronal contribution, like- 
wise there is evidence for high density plasma at cool tem- 
peratures. SU Aur is the by far brightest and with tem- 
peratures of at least up to 50 MK hottest X-ray source in 
our sample. A powerful and active corona is the essential 
contributor to its X-ray spectrum and frequent flaring is 
observed. Unfortunately the very strong absorption pre- 
vents definite conclusions about the cooler plasma com- 
ponent. TW Hya is the other hand strongly dominated by 
accretion, but an additional coronal component is clearly 
detected. It is the prototype and still the outstanding and 
only example of an accretion dominated star. Most of its 
plasma is at cool temperatures typical for accretion spots 
and it exhibits by far the lowest O vn f/i ratio. Also the 
highest neon to oxygen ratio, which may be interpreted via 
depletion of grain forming element, is found for TW Hya. 
This is plausible since TW Hya is the oldest and most 
evolved object in our sample. In the other sample CTTS 
neon is also enhanced, but the abundance pattern is more 
reminiscent of the inverse FIP effect found in active stars. 

Considering the global picture of stellar evolution to- 
wards the main sequence, we find that magnetic processes 
play a major role in high energy phenomena in all our 
sample CTTS. While the specific stellar evolution and its 
timescale might depend on e.g. spectral type, the genera- 
tion of X-rays via accretion and coronal activity is appar- 
ently a common feature of CTTS in general. 
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